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Abstract. - We consider ferromagnetic particles embedded into a superconductor and study 
the screening of their magnetic moments by the spins of the Cooper pairs in the superconductor. 
It is shown that a magnetic moment opposite to the one of the ferromagnetic particle is induced 
in the superconductor. In the case of a small itinerant ferromagnet grain and low temperatures 
the full screening of the magnetic moment takes place, i.e the absolute value of the total 
magnetic moment induced in the superconductor is equal to the one of the ferromagnetic 
particle. In type II superconductors the proposed screening by spins of the conduction electrons 
can be much stronger than the conventional screening by Meissner currents. 



The phenomenon of screening is very common in physics. The best known example is 
the screening of an electric charge in metals due to a redistribution of free electrons in space. 
This charge screening is very strong and the length characterizing an exponential decay of 
the electric field (the Thomas- Fermi length \tf) is in most metals of the order of the Fermi 
wave length, i.e. of the order of the interatomic spacing. The length Atf does not depend 
on whether the metal is in the normal or in the superconducting state. 

Another famous example is the screening of a magnetic field or a magnetic moment by 
superconducting currents in a superconductor (Meissner effect) [1]. Due to this effect the 
magnetic field decays over the length (the penetration depth) and vanishes in a bulk 
superconductor. The same length characterizes the decay of the magnetic field created by a 
ferromagnetic {F) grain embedded in a superconductor. In contrast to the screening of the 
electric charge in normal metals the screening of the magnetic field and magnetic moments 
in superconductors is weaker and the penetration depth Al can be of the order of hundreds 
interatomic distances or larger. The screening of the magnetic moment is a phenomenon 
specific for a superconductor and, in contrast to the charge screening, is very small in a normal 
metal. Although the stray field created by a ferromagnetic grain embedded in a nonmagnetic 
{N) metal may induce a negative local magnetization in some regions of the normal metal due 
to the Pauli paramagnetism, the susceptibility is rather small {ji^gV ~ 10~^, /is is the Bohr 
magneton and v the density of states at the Fermi level) and the screening can be neglected. 
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For certain geometries of a superconducting sample (films, wires) the penetration length Al 
can exceed the transversal size and the screening of the magnetic moment of a ferromagnetic 
particle due to the Meissner currents does not play an essential role. It is usually believed 
that in such a situation the total magnetic moment is just the magnetic moment of the 
ferromagnetic particle and no additional magnetization is induced by the electrons of the 
superconductor. This common wisdom is quite natural because, at first glance, a possible 
contribution into the screening of the electron spins in the superconductor is even smaller 
than in the normal metal. In conventional superconductors the total spin of a Cooper pair 
is equal to zero and the polarization of the conduction electrons is even smaller than in the 
normal metal. Spin-orbit interactions may lead to a finite magnetic susceptibility of the 
superconductor [2], but it is positive and anyway smaller than in the normal state. 

In this Letter we suggest a new mechanism for the screening of the magnetic moment 
of a ferromagnetic particle embedded in a superconductor by spins of the superconducting 
electrons. This effect is large and the magnetic moment of the grain can be completely 
screened. The characteristic length of the screening is of the order of the size of the Cooper 
pair = Ds (we consider the "dirty" limit) and can be much smaller than the 
penetration depth Al in type II superconductors. If the size of the superconductor in the 
transverse direction is smaller than the penetration length Al, the mechanism we propose is 
the only one leading to the screening of the magnetic moment of the ferromagnetic grain. 

Although this effect seems very surprising and has been overlooked in all previous inves- 
tigations of superconductivity, its origin can be understood without any calculations. This 
additional screening arises due to the exchange interaction between the spins of the conduction 
electrons and the magnetic moment of the F grain, and the possibility for the superconducting 
condensate to penetrate the grain. 

If the size of the grain is much smaller than the size of the Cooper pair ^5 the probability 
that both electrons of a Cooper pair are located in the ferromagnetic grain is small. Therefore 
one can assume that only one electron of the Cooper pair spends some time in the grain. Then, 
the exchange interaction enforces the spin of this electron to be parallel to the magnetization in 
the grain. This leaves no choice for the second electron in the Cooper pair but to be antiparallel 
to the magnetization of the grain. In this way an additional magnetization antiparallel to the 
one in the grain is induced in the superconductor. Of course, if the transparency of the 
S/F interface is small the induced magnetization is weak. However, if the transparency is 
high enough the induced magnetization is large and the magnetization of the grain can be 
completely screened. 

Below we support this qualitative discussion by an explicit calculation using quasiclassi- 
cal Green's functions and a proper Usadel type equation including the exchange field. Al- 
though the calculational scheme is quite standard, it contains a new important ingredient: 
the superconducting condensate in the F region near the surface consists of a singlet and 
a triplet components [3]. Due to the proximity effect the triplet component penetrates the 
superconductor over the length ^5 and, as it carries spin, polarizes the superconductor. The 
triplet component was not considered previously for this type of problems. We consider a 
ferromagnetic spherical grain of the radius a embedded in a superconductor S (Fig^). The 
superconductor is assumed to be a conventional s-wave superconductor with singlet pairing. 
The Hamiltonian describing the system can be written as 



where Hq is the one-particle electron energy including interaction with impurities, J is the 
exchange field which is nonzero only inside the F grain, and A is the superconducting order 




(1) 
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Fig. 1 - Schematic representation of a magnetic grain embedded in a superconductor. The arrows 
show the direction of magnetization. In the superconductor a negative magnetization is induced over 
distances of the order of ^s- 



parameter which vanishes inside F. The matrix 173 is the PauH matrix in the spin-space and 
we have assmned that the magnetization of the grain is in the z direction. The notation p 
and s means inversion of momentum and spin respectively. 

In order to calculate physical quantities of interest it is convenient to use quasiclassical 
Green's functions 5 and write the proper quasiclassical equations for them. We consider the 
"dirty" limit, which means that the inverse momentum relaxation time r^^ is much larger 
than the exchange energy J and than the critical temperature of the superconductor. The 
functions g are 4x4 matrices in the particle-hole and spin space and obey the Usadel equation 
which in the presence of the exchange interaction has the form 

^V(5V5)-c^[f3<To,g]+zJ[f3(73,.g] = -z[A,g] . (2) 

In the S region D — Ds, J = 0, A — Aif2(73 (the phase of A is chosen to be zero). In the F 
grain D = Dp and A = 0. Eq. |2l should be complemented by boundary conditions at the 
S/F interface [4] 

IF {9n\7g)p = 75 (snVg)^ ; 7f {gn\/g)p = - [53, 5^] , (3) 

where 'ys,F = o's.pRt, <^s,f are the conductivities of the F and S layers, Rh is the resistance 
per unit area, and n is the unit vector normal to the S/F interface. Solving Eqs. l(21 13|) for 
the Green functions g, one can obtain the induced magnetization SM 

cj— +00 

6M = ^6^{< c|,^Cp| - c^^Cpi >) = -^iiTTvT ^ Tr{a3g)/2 (4) 
p UJ——00 

where /i is an effective Bohr magneton, ly is the density of states (DOS) and the sum is taken 
over the Matsubara frequencies oj — ■nT{2n +1). 

The functions g can be represented in the form g = if^j -f Tsg, where the condensate 
function / and the normal function g are matrices in the spin space [3]. For simplicity we 
consider the case of a small grain when the condition a < = yDp/J , fulfilled (this is 
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possible provided the exchange energy J is much smaller than the Fermi energy ep)- In this 
limit the solution can be found averaging Eq. Q over the grain volume and we obtain for the 
Green functions 

gF± = C^±/(lj±, fF± ^ ±ebFfBCs/Cuj± (5) 

Here gp± and fp± are the diagonal elements of the matrices g and /, Lij±=uj + ebrgscs T ^J, 



Clj± = -\/w± - [ibpfBCsY, 9BCS = i{uj/A)fBcs w/Vt^^ + A2, EbF = 3DF/{2jFa). When 
obtaining Eq.©, we assumed that the functions gs, fs are close to their BCS values gscs-, 
Jbcs, i-e- the corrections Sgs, Sfs are small. Under this assumption Eq.(|3) for the functions 
6gs, Sfs can be linearized. For example, the linearized Usadel equation for Sgss = Tr{asg)/2 
has the form 

V'(5gs3 - 45gsi = (6) 

The function Sgsz determines the excess magnetization. Solving Eq. with the boundary 
condition, Eq. (j2Jl, we obtain 

IbCS ^ g^^s(^~*i) 

Sgs-i = {gBcsfFo - fBcsgFs) -r- , (7) 

75 1 + Ksa r 



where k| = 2\/uj^ + ^ / Ds, Jfo = {fF+ — fF-)/2 describes the triplet component mentioned 
above. One can see that the correction Sgss is small if the parameter a/js is small. Note also 
that the function Sgs3, which according to Eq.0) determines the induced magnetization, is 
proportional to the condensate function fpQ. The latter function is the triplet component with 
zero projection of the magnetic moment on the z axis. Using Eq.Q) one can easily calculate 
the total magnetic moment of the S region 



- + OC 



U — — OQ 



Ms^ -iniysTfi > / d-'rdgss ■ (8) 



It is not difficult to see that the magnetic moment Ms has the sign opposite to J (the 
magnetic moment of the ferromagnetic particle Mpo is proportional to J) , which means that 
the spins of the conduction electrons screen (at least partially) the magnetic moment Mpo of 
the ferromagnetic grain. 

Let us compare the total magnetic moment induced in the superconductor, Eq.®, with 
the magnetic moment of the ferromagnetic grain in the normal state Aipo = {Ana^ /3)MpQ. 
The sum in Eq.® can be computed numerically in a general case. For simplicity we consider 
here a limiting case that can be realized experimentally. 

We assume that the transmission coefhcient through the S/F interface is not small and 
the condition A << J < [Dp /a?) is fulfilled. In this case the expression for fpQ is drastically 
simplified. To estimate the energy Dp/ a? we assume that the mean free path is of the order 
of a. For a = 30A and vp = lO^cm/ sec we get Dp/a^ « lOOO-fC;. This condition is fulfilled 
for ferromagnets with the exchange energy of the order of several hundreds K. In order to 
relate MpQ to J, one has to make a certain assumption about the nature of the ferromagnet. 
If the magnetic moment Af^o is induced mainly by free electrons (an itinerant ferromagnet), 
one gets MpQ — fivpj. Then we obtain for low temperatures 

Ms/Mpo = -1 . (9) 

Eq. lO describes a remarkable phenomenon: at sufficiently low temperatures and in the limit 
of a small grain {a < ^p) the magnetic moment of the latter is screened over distances of the 
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order of ^5. This screening is complete if the magnetization of the F particle is due to the 
free electrons (itinerant ferromagnet). It can be easily shown that at arbitrary temperatures 
this ratio is equal to —(1 — n„(T)/ne), where rie and n„(T) are the density of total number 
of electrons and " normal" electrons defined in Ref. [1] . 

The compensation of the magnetization of an itinerant ferromagnet by the Cooper pairs 
is to some extent consistent with the result obtained by Rusinov and Gor'kov some decades 
ago [5]. They studied the properties of a superconductor with paramagnetic impurities that 
were assumed to be ferromagnetically ordered. The free electrons interact with the magnetic 
impurities via the exchange interaction. Their approach (averaging over impurity positions) 
reduces the problem to finding the magnetization of a superconductor with an effective ex- 
change interaction uniformly distributed in space. It was demonstrated in Ref. [5] that the 
total itinerant magnetization of the system was zero in the limit of low temperatures and not 
too large exchange energy. 

Clearly, the screening obtained here is due to the appearance of the triplet component. 
So, if the latter is suppressed by other mechanisms the effect will be reduced. The spin-orbit 
interaction (SOI) and orbital effects (Meissner currents) are such mechanisms. The other way 
of thinking of this reduction is that the spins of the electrons of the Cooper pairs are not 
necessarily antiparallel in the presence of e.g. SOI and are not as efficient in inducing the 
magnetization. 

In Ref. [3] we studied the effect of the SOI on the triplet component penetration into the 
ferromagnet. Here we can simply use these results to analyze the effect of the SOI on the 
penetration of the induced magnetization into the superconductor. 

If the SOI is taken into account, an additional term of the form {i/Tg.o.) [ST^gT^Sjg^ 
appears in the Usadel equation, where Tj o. is the spin-orbit scattering time, S — (cri, a2, 173X3) 
(see [6] and [3] ) . Due to this additional term the quantity k| in the linearized Usadel equation 
is replaced by 

4 = 4 + 4o , (10) 

where = 8Ds/tso- Therefore the length of the penetration of fsa and of Ms into the S 
region decreases if k| ~ < n^^. In principle, one can measure the spatial distribution of 
the magnetic moment in the S region (see e.g. [7] ) and get an information about the SOI 
in the superconductors. With the help of Eq. ifTUI) . this would be an alternative method 
to measure the strength of the SOI in superconductors, complementary to the measurement 
of the Knight shift [8]. (The latter is based on the result of Ref. [2] that the Knight shift 
observed in superconductors is due to the SOI). Let us notice, that in the presence of SOI the 
linearized equation for the function fss that determines a correction to the energy gap due to 
the proximity effect remains unchanged. 

The orbital effects (the Meissner currents) also change the characteristic length of the 
penetration of the induced magnetization Ms- Contrary to the case of bulk conventional 
superconductors, the Meissner currents in S/F structures arise spontaneously even in the 
absence of an external magnetic field Hext- This happens because an internal magnetic field 
is induced by the ferromagnet. These spontaneous currents were studied, e.g., in Refs. [9,10]. 
In Ref. [9] the spatial dependence of the Meissner current in the F film was calculated, whereas 
the Meissner currents in the both S and F films were computed in Ref. [10]. It is of interest 
to know not only the spatial dependence and the magnitude of the Meissner currents, but also 
their effect on the penetration length of the TC and on the induced magnetization. Here we 
study the infiuence of the Meissner currents on the penetration length of Ms- For simplicity 
we consider first a planar geometry, which allows us to get a simple solution. For a spherical 
particle we estimate the effect by order of magnitude. 
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Let US consider a bilayer S/ F structure with the thicknesses ds^p- We assume again that 
the thickness of the F layer is small {dp << £,f ) and the thickness of the S layer obeys 
the condition: £,s << ds << Xl- In this case the solution for gp± and fp± (Eq.©) 
remains unchanged provided one replaces a/3 with dp . The solution for Sgs3 has the form 
Sgs3 = ii'isgBCsfBCs/ls)fposiip{~Ksx) and one can easily calculate the spatial distribution 
of the magnetization in the system Ms,p{x). Note that again in the case of an itinerant 
ferromagnet and low temperatures the total magnetic moment in S is equal to (— Mp^dp). 
We display schematically the spatial dependence of the magnetization in Fig. 2 alongside with 
the spatial dependence of the vector potential A{x) which is given by 

A[x) = + 47r ( dx'M{x'). (11) 

with a constant determined from the condition that the total current through the system 
is zero (no external magnetic field). We have assumed that the magnetization lies in-plane. 
The supercurrent density is expressed through A as 

3{x) = 27rTaA(a:)/(e<^o) ^3 (12) 

where (^o is the magnetic flux quantum. In the limit J « the condensate functions /I 
are nearly the same in the S and F layers. Taking this into account and calculating the total 
current /, we get from the condition / = for Aq 

asds + apdp/2^as—J2 /i'^sds + apdp) . (13) 

where A„i — AirMpodp . At low temperatures the third term is approximately equal to 
0.53as\/ Ds/2A. The spatial dependence of the vector potential is shown in Fig. |21 Eq.© 
for 6gs3 does not change in the presence of the Meissner current provided one replaces 
k| with Kg + Pq, where po — Ao/(j>o. It is clear that the orbital effects are negligible 
if Ao/(t)Q = Am/fj^o « Cs^- For example, Am/cfio ^ S.lO^cm^^ for AirMp ~ IkOe and 
dp ^ 50A, and therefore should be smaller than 2/im. 

It is not difficult to estimate po in the case of a spherical particle. The quantity pq in this 
case reaches the maximum at r = a and is of the order of {'ii:Mp / 4>o)a- As we have noted 
above, the magnetic moment in F is compensated by that in the superconductor. This means 
in particular that the total internal magnetic flux in the S/F/S Josephson junction may be 
zero even if the F layer is a single domain (see experiments of Ref. [11]). 

In summary, we have shown that the magnetic moment MpQ of a ferromagnetic particle 
embedded into a superconductor is screened by spins of the Cooper pairs. An induced magnetic 
moment Ms aligned in the opposite to AIpQ direction arises in the superconductor. 

The induced magnetic moment in the superconductor can be observed experimentally. 
In a recent experiment [7] the spatial electron spin polarization was determined by means 
of muon spin rotation. Such a method may be used in order to determine Ms- Another 
possible method is the measurement of the Knight shift in superconductors that should be 
dependent on the magnetization Ms -One can also determine the total magnetic moment of a 
S/F structure performing magnetic resonant measurements as in Ref. [12]. 

At last, one can determine Mg measuring the total magnetic moment of a superconductor 
with embedded ferromagnetic particles similarly to the work [13], where an enhancement of 
the magnetic moment of ferromagnetic particles embedded into a nonmagnetic matrix has 
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Fig. 2 - Spatial dependence of the vector potential A{x) or the supercurrent density. For low enough 
temperatures the values of A at both outer surfaces {x = —ds and x — dp) is nearly the same. The 
dot-dashed line shows the spatial dependence of the local magnetization in arbitrary units 

been observed. Although the enhancement in the case of a normal metal matrix awaits its 
explanation, we suggest to measure the magnetization of the system of ferromagnetic particles 
embedded in a superconducting matrix. When using a superconducting metal as the matrix, 
a reduction of the magnetic moment, instead of the enhancement, should be observed below 
the superconducting temperature Tc- 
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